Quantum Down Conversion and Multipartite Entanglement via a Mesoscopic SQUID 

Ring 
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In this paper we study, by analogy with quantum optics, the SQUID ring mediated quantum 
mechanical interaction of an input electromagnetic held oscillator mode with two or more output 
oscillator modes at sub-integers of the input frequency. We show that through the non-linearity of 
the SQUID ring multi-photon down conversion can take place between the input and output modes 
with the resultant output photons being created in an entangled state. We also demonstrate that 
the degree of this entanglement can be adjusted by means of a static magnetic flux which controls 
the strength of the interaction between these modes via the SQUID ring. 



I. INTRODUCTION 

Over the last five years there has been a rapidly grow- 
ing interest in quantum technologies, particularly for the 
fields of quantum information processing and quantum 
computing^. This has been aided both by the creation 



of key algorithms^ and by improvements in experimen- 
tal technique. A range of approaches, involving different 
physical disciplines, have been adopted in the pursuit of 
these technologies but at the basis of each is the concept 
of using entangled states of quantum objects^. The cre- 
ation of entangled states, and the exploitation of their 
properties, are therefore seen as fundamental to success- 
ful progress in these technologies. In particular, it is 
clear that there exists a need to develop systems or de- 
vices through which the degree of entanglement between 
quantum objects can be actively controlled. In this paper 
we show how the intrinsic non-perturbativc behaviour of 
a quantum SQUID ring (here, a thick superconducting 
ring enclosing a Josephson weak link device) can be used 
to achieve this goal, taking as our specific example en- 
tanglement between photon states. 

Interest in the use of SQUID rings (and other weak 
link circuits) as quantum logic devices has been stim- 
ulated by recent experiments to probe superposition of 
states in these systems£*£*i2*2*iS. Furthermore, other re- 
cent experiments have been performed that show coher- 
ent dynamics of superconducting flux qubits coupled to 
individual harmonic oscillators^. Building on this back- 
ground, and previous theoretical work by us on the cou- 
pling of SQUID rings to electromagnetic (em) fields^, 
we consider the process of photon down conversion via 
a mesoscopic SQUID ring. Specifically, by analogy with 
the field of quantum optics, we demonstrate that the ring 
may be viewed as a non-linear medium which can down 
convert photons from one frequency (input field oscillator 
mode) to another lower frequency and generate entan- 
gled photons (output field oscillator modes), as is often 
performed in quantum optics experiments^. An analy- 



sis of multiple output modes is important since it allows 
us to understand the way in which several quantum cir- 
cuits may be entangled - an important issue in quantum 
technologies - for example, as entanglement registers in 
quantum computing. We also note that by coupling in a 
larger number of field modes we change the position (in 
$ x ) of the splittings in the spectrum of the Hamiltonian 
for the composite system. We consider, therefore, that 
it is by no means trivial to assume that we can simply 
extend the results of our previous worki^iil. As we have 
shown by computation, in a three mode system (two em 
field modes plus one SQUID mode) entanglement can 
be adjusted through the control of a static bias flux Q x 
applied to the ring. From an experimental viewpoint, in 
this system down conversion processes are relatively easy 
to observed at least for the situation where microwaves 
(input electromagnetic field as in the range of a few to 
10GHz) are down converted to radio frequencies (rf out- 
put mode as 1 to 20MHz)i£. We note that microwave 
photons are difficult to transport, which could be prob- 
lematic for some applications in quantum communica- 
tion. By implication, but not calculated in this work, 
it should also prove possible to observe up conversion 
processes involving significant frequency ratios between 
input and output. However, the ultra low noise elec- 
tronics required for the output stage in an experiment of 
this type (say, for inputs at rf and outputs at microwave 
frequencies) is much more difficult to engineer than at 
the much lower radio frequencies. Having observed down 
conversion in the experiments at Susse:si£, we tried to 
emulate this behaviour theoretically. However, such large 
frequency differences (a few GHz input, a THz SQUID 
ring characteristic oscillator frequency and a few MHz 
output) was found to be computationally intractable at 
this time. Bearing this in mind, and considering the in- 
terest currently shown in quantum computing, we chose 
to demonstrate in this paper that SQUID rings can act 
as devices which mediate large ratio frequency down con- 
version that is also computationally tractable. Moreover, 
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for multiple output modes we demonstrate that the down 
converted photons are entangled together. In our com- 
putations we have modelled the field oscillator modes 
as equivalent resonant LC circuits. Experimentally, this 
could be realised by utilising real LC circuits constructed 
from superconducting material to minimise decoherence. 
Alternatively we could use a transmission line comprising 
a series of coupled superconducting LC rcsonator o 15 ! 16 . 

From an experimental perspective, it is interesting to 
note that there has already been some progress made in 
realising simple quantum circuits systems relevant to the 
work presented in this paper. In this regard there are 
now several examples of coupled superconducting qubits 
(here the SQUID ring is, to all intents and purposes, a 
flux qubit). These include fixed coupling flux qubits^, 
current biased Josephson-junction qubits that are either 
coupled capacitivelj^ or via a connecting superconduct- 
ing loopiS. Recently, there has even been an example 
of tri-partite quantum entanglement in macroscopic su- 
perconducting circuits^. Moreover, an application that 
relates very strongly to the circuit systems presented in 
this work utilised a SQUID-ring to mimic a simple har- 
monic oscillator. This device was then coupled, quantum 
coherently, to a superconducting flux qubi^ii. However, 
although the physics of all these systems is similar, in 
this work we have chosen to consider a superconduct- 
ing circuit coupled to quantum electrodynamic degrees 
of freedom, i.e. photons in a cavity. Recently it has been 
shown that a superconducting qubit can be quantum co- 
herently coupled to an on-chip cavityi^. Furthermore, we 
note that superconducting devices, which can be coupled 
to these and more exotic systems, such as GHz nanomc- 
chanical resonators^, are not only easy to fabricate, but 
are adaptable to a wide variety of applications. 

In this work we first investigate, from a full quantum 
mechanical viewpoint, the simplest down conversion and 
entanglement process, involving one input photon at a 
given frequency and two output photons at half this in- 
put frequency. To the limits of the computational power 
available to us, we extend our theoretical investigation 
of such down conversion processes to that of generating 
four output photons, each at a quarter of the frequency 
of the input photon (i.e. a four photon down conver- 
sion). Within the Adami-Cerf criterion22i2&24 we show 
that entanglement exists between any one of these output 
oscillator modes and the rest of the system (i.e. SQUID 
ring and input oscillator mode) . This implies that all the 
output modes will be entangled with each other, as well 
as to the rest of the coupled system. While this multipar- 
tite entanglement may not be as strong as some bipartite 
entanglements, such as Bell States, it is a key component 
in certain requirements in quantum computation, such as 
entanglement registers, where a large number of qubits 
are entangled together to generate a total entangled state 
for the process of computation. It is this idea that is re- 
sponsible for the (as yet in principle) power of quantum 
computing over modern day, classical computers. Also, 
the question of characterising entanglement in multipar- 



tite quantum systems is still very much an open one. 
Whilst each measure agrees as to whether the state vec- 
tor of a given system is separable, or maximally entan- 
gled, they do not, in general, provide an ordering that 
can be used to quantify entanglement of more general 
states. In this paper we have chosen to use an entangle- 
ment measure that gives us an idea of the distribution of 
the quantum correlations that exist in our example mul- 
ticomponcnt systems. However, as we have not produced 
maximally entangled states it is not possible to quantify 
absolutely the level of entanglement present. We note 
that for a simpler system, controlled via a dynamic flux 
biasS^, we have shown that it is possible to create a max- 
imally entangled state between a SQUID ring and an 
em-field mode. It may be that similar methods could be 
used in these composite systems to produce such states 
of maximal entanglement. Unfortunately computations 
with time dependent Hamiltonians are far harder to solve 
than the static cases presented in this paper, and are be- 
yond the limit of current computational power. 

For the example of the two photon process we consider 
the effect of dissipation on the coherent evolution of the 
coupled (entangled) system. By extension, we show by 
computation that much higher order (i.e. more than our 
computed four photon case) down conversions may prove 
possible experimentally To this effect we provide a fur- 
ther example of down conversion by a factor of twenty 
between the input em oscillator mode and an output 
mode. However, both due to limitations in computer 
power, and a still insufficient understanding of entangle- 
ment in many particle (multipartite) quantum systems, 
we have not proceeded any further in characterising this 
high order example. 



II. THEORETICAL MODEL 

In our theoretical treatment of the interaction of 
quantised em oscillator modes (photons) with a meso- 
scopic quantum mechanical SQUID ring wc model the 
em field modes as parallel LC resonators, as a sim- 
ple and convenient way of parameterizing the (physi- 
cal) cavity modes coupled to the SQUID ring. The 
scheme we have adopted is shown diagrammatically in 
figure n for the simpler, two photon, down conver- 
sion process. In this example the input photon fre- 
quency is w ei /27r while the two (identical) output fre- 
quencies are w e2 /2i: and w e3 /27r such that uj e . 3 = uj e3 = 
0.5o;ei- Here, the SQUID ring characteristic oscilla- 
tor frequency is lj s /2tt = l/2it\/ A S C S for a ring induc- 
tance and a weak link capacitance of A s and C s , respec- 
tively. In our calculations we treat this system fully quan- 
tum mechanically by assuming an operating temperature 
T such that hw s (= h/^A s C s ) , hw ei (= hj */ C ei L ei ) and 
hoj e2 3 (= %/ \JC e23 L e2 3 ) are all much greater than fcgT 
for field mode inductances and capacitances L ei and C ei 

(where, generally, i — 1,2,3 ). In this model system 

the coupling between the SQUID ring and the oscillator 
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FIG. 1: Schematic of a four mode system comprising three 
electromagnetic field mode inductively coupled to a meso- 
scopic SQUID ring. 



modes is taken to be inductive. As can be seen in fig- 
ure Q] we have chosen to couple these circuits so that all 
interactions between the field modes are performed via 
the SQUID ring. In this four mode system (i.e. three em 
modes plus one SQUID mode) the static flux $ x applied 
to the ring can be used to control the level of interaction 
between the input em mode and an output em mode of 
the system. In a specific example in this paper, as an 
extension of previous worki^, we demonstrate that the 
behaviour of coupled SQUID ring-em mode systems is, 
in general, very strongly dependent on this control flux. 
Although this aspect of our investigations is an exten- 
sion of the results of our previous papersi2ii4, it is im- 
portant to show that this dependence remains valid for 
the kind of multicomponent circuits that would appear to 
be required in the quantum technologies currently being 
pursued 1 . 2 . 5 . 6 ! 7 ! 8 ! 9 ! 10 . 11 ! 24 . 

The Hamiltonian for the coupled system of figure ^ is 
made up of the uncoupled Hamiltonians for each mode 
together with a set of interaction terms. For each em field 
mode (both input and output) the Hamiltonian takes the 
formic 



H, 



<I>: 



2C K 



2A f 



(1) 



where, again, i = 1,2,3 , while for the SQUID ring we 

use the standard, time independent, lumped component 
circuit HamiltonianSLSl (the wavelengths of the photon 
modes being very much larger than the dimensions of the 
mesoscopic SQUID ring circuit). Specifically 
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(2) 



In this description <& s (the magnetic flux threading the 
SQUID ring) and Q s (the electric displacement flux be- 
tween the electrodes of the weak link) are the conju- 
gate variables for the ring, hv/2 is the matrix element 



for Josephson pair tunnelling through the weak link (of 
critical current I c = 2ev) and $o = h/2e is the su- 
perconducting flux quantum. From the perspective of 
the work presented in this paper, the Josephson cosine 
term in the SQUID ring Hamiltonian J3J allows for non- 
perturbative behaviour to all orders in the coupled sys- 
tem of figure n I n turn, this means that in interactions 
between a quantum SQUID ring and external em fields, 
non-perturbative, multiphoton absorption/emission pro- 
cesses tend to dominate within an exchange regioni£. 

To illustrate the form of the solutions to the time in- 
dependent Schrodinger equation for a SQUID ring, using 
the Hamiltonian we show in figure Ufa) the first four 
cigenencrgies (n = 0, 1, 2, 3) as a function of external flux 
<E> X for a typical mesoscopic quantum SQUID ring with 
parameters^ C s = 5x l(r 15 F, A s = 3 x l(r 10 H (giv- 
ing a characteristic SQUID ring oscillator frequency of 
130GHz, this is much smaller than the superconducting 
energy gap in, for example, Niobium, that is often used to 
fabricate SQUID circuits) and hv = 0.035$q/A s equiv- 
alent to a Josephson frequency of 700GHz, which is not 
unrealistic when compared to those used by Friedman 
et al 6 . (2THz) and Orlando et al2& (2.12THz). These fre- 
quencies should be sufficient to avoid thermal excitations 
for any experiments performed at 40mK, a temperature 
easily obtainable in modern dilution refrigerators. If wc 
assume that the Josephson weak link in the SQUID ring 
is of the tunnel junction type with, typically, an oxide 
tunnel barrier of dielectric constant 10, then the junc- 
tion dimensions will be close to 0.2/im square, which can 
be fabricated using current lithographic and thin film 
techniques. Furthermore, with hv = 0.035<I>o/A, and 
these dimensions, the corresponding Josephson critical 
current density is w 4kA/cm 2 which appears suitable 
for SQUID rings operating in the quantum regime. In 
addition to the ring eigencnergies, we also show in fig- 
ure a) the harmonic oscillator field mode cigenencrgies 
of ei (huj ei = Hu) s /2) (dotted) and e2,3 (huj e23 = Hw s /4) 
(dashed). In figure Gib) we show the eigenenergies of the 
total coupled system. We note that at the level of capac- 
itance chosen for the SQUID ring the energy splittings 
where the field mode energies and the ring energies cross 
are very small. In figure Efb) almost all these splittings 
are unresolved on this scale but at sufficient resolution 
all would be discernible. 

Provided that we disallow any direct interaction be- 
tween the field mode oscillators (this being our assump- 
tion and a situation which is easy to establish experimen- 
tally), the coupling terms are given by 



where the fi eiS are the fractions of magnetic flux coupled 
inductively between any one of the field mode oscillators 
and the SQUID ring. The total system Hamiltonian H t 
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FIG. 2: (a) the eigenenergy spectrum of an uncoupled mesoscopic SQUID ring at T = OK as a function of external bias flux 
(solid lines) together with the harmonic oscillator spectra for the input mode (dotted) and the two output modes (dashed) 
(b) the eigenenergy spectrum of the coupled system (figure ^ for the four modes coupled together inductively with the circuit 
parameters given in the text. 



for the three em mode system then reads 

3 

H t = H s +J2(H ei -H eiS ) (3) 
i=i 

We can make a transformation using the unitary oper- 
ator T = exp (—i& x Q s /h) to translate the SQUID Hamil- 
tonian H s into a more convenient form. The Hamiltonian 



(0) can then be written asi£ 

where the "I^-periodic dependence of the (transformed) 
ring Hamiltonian on $> x is explicit. This unitary transfor- 
mation also modifies the ring-field mode interactions to 
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the form H' es = !j^-<& s <& ei so that, following on from 
the transformed system Hamiltonian (again $o-pcriodic 
in $ x ) can the be written as 



Hi = H s 



Adopting this Hamiltonian we can solve the cigcnprob- 
lem H[ |£„) = S„ |£„) and use the eigenvectors |£„) to 
construct the evolution operator U(t) for the systemic 
via, 



(tn\ 



(5) 



The time averaged energy expectation values {(Hi)) for 
the ring and the field modes can then be calculated from 
the expression 



((Hi)) = lim / Tr [ Pi (t) H t ] dt 



(0) 



where the pi (t) are the reduced density operators for 
the sub-components i = ex,e2,e3 and s for the coupled 
system^ which are obtained from the density operator 
p(t) = U(t)p(0) (for an initial density operator p(0)) by 
using pi (t) = Tr^i {p (t)}. 



A. Four mode (two photon) down conversion and 
entanglement 

1. Without dissipation 

Referring again to the system of figure ^ an d assum- 
ing that it starts in an initial state |^/>(0)}, the evo- 
lution operator Q can then be used to generate the 
system state vector \ip (t)) = U(t) \ijj (0)) at any subse- 
quent time t (from which p(t) can be found via p(t) = 
|?/> (t)) (ip (t)\). As an illustration we now assume that 
the system, SQUID ring + oscillator modes, starts in the 
state \ip(0)) = |l) ei ® \ a ) s ® |0) ea <g> |0> ea , where a is 
the ground state of the ring and the \n) are the photon 
number states for the various field modes of the system. 

Starting in this initial state we show in figure yfa) 
the time averaged energy expectation values of the em 
field modes (ei - the input and e2.3 - the two outputs) of 
the system as a function of 4> x , where the two outputs 
(e2 and e^) superimpose exactly, as is to be expected by 
symmetry. As is apparent, due to the interaction of the 
SQUID ring with e\, and one of the other two (output) 
modes e2 and e^, sharp transition regions (or energy ex- 
change regions^) develop in the time averaged energies, 
discussed by us in a previous publication^. These re- 
gions of energy exchange between the ring and the field 
modes occur over very narrow ranges in <& x , which is 
clear from our example system in figure |3fa) and (b). As 
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FIG. 3: (a) the time averaged energy expectation values for 
the electromagnetic input and output modes of the coupled 
system of figure [5] as a function of external bias flux (b) an 
zoomed in enlargement of the central section of the results 
in (a) as a function of bias flux showing a series of energy 
exchange regions; in the calculations presented subsequently 
the two bias points A and B have been selected, the former 
inside an exchange region, the latter outside. 



we shall demonstrate, the non-perturbative nature of the 
ring in these exchange regions is sufficient to generate 
photon down conversion and entanglement between the 
field modes. In an earlier calculation^^ we showed that 
ring-field mode interaction strength is a maximum at the 
centre of the exchange region and very small at its edge. 
For the calculations presented in this section we have 
chosen to bias the coupled system at the centre of the 
exchange region, denoted A in figure |3^b). The energy 
exchange region is seen as the energy splitting (divided 
crossing) circled in figure |2Ib). 

In figurcs^Ja) and (b) we show how the number expec- 
tation values (n ei ), (n e2 ) and (n es ) of the three em field 
modes vary with the dimcnsionless time (r = to s t) for the 
bias points A and B in figure|2Ib), where Q x = 0.49I83$o 
(at A in figure H») and $ x = 0.487$ (at B in fig- 
ure^Ib)). Here, we note first that for bias point A there 
is down conversion of energy from field mode ei to one 
or other of the two half frequency (w ei /2) modes e2,3 
and second that (n e . 2 ) and (n e3 ) versus uo s t are iden- 
tical due to the symmetric nature of their coupling to 
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grey in figure^Ja). As expected, with e\ starting in state 
|0) ei there was no down conversion. However, we note 
that there is some occupancy in the output modes due 
to the ring-mode coupling energy available in the sys- 
tem. It is, of course, possible that higher harmonics of 
the output field oscillators could be excited, which may 
not be desirable. In order to investigate such excitations 
we have included an additional output mode, at the same 
frequency as the input mode, and recomputed the results 
of figure Ufa). This is shown in figure 0Jc) Although it is 
apparent that the level of down conversion is reduced it 
is, however, still significant. It would appear, therefore, 
that the presence of unwanted higher harmonics in the 
output modes is a factor which might need to be reduced. 
Nevertheless, it does not negate the computations where 
such processes are not considered. 

Concomitant with the two photon down conversion 
ei — ► 62,3 there exists the possibility that the photons 
in the lower frequency modes e2 and e% are strongly cor- 
related or even entangled. In order to characterise the 
correlations within this system we will employ two en- 
tropic quantities. The first of these yields information 
about the correlations that exist between any two com- 
ponents of a systemi2iSL2ii22iS and is defined by 



I(A,B) = S a + S b -Sai 



(7) 



DimciKlonldsa Time 



where S m is the von Neumann entropy for each subsys- 
tem (m)2i. We note that for a bipartite closed system 
this measure is often employed to characterise the en- 
tanglement present. However, complications arise when 
this quantity is used in multipartite or open systems (for 
more information see (34|). Nevertheless, it is still useful 
for gaining an appreciation of the correlations that are 
present between the components of the system in which 
we are interested. The second entropic quantity we shall 
use is the entanglement measure introduced by Adami 
and CerS£, which for a subsystem A of the total system 
T is given by 



FIG. 4: (a) number expectation values for the electromag- 
netic input and output modes of the coupled system of fig- 
ure |2| as a function of dimensionless time, at external flux 
point A in figure starting with a photon occupancy of 1 in 
the input mode; also shown in grey is the output mode num- 
ber expectation values for the zero photon occupancy in the 
input mode (b) the same results as (a) above calculated at 
the external bias flux point B in figure (c) a repeat of the 
calculation of (a) with the inclusion of a third output mode 
e4, with its number occupancy shown in grey, at the same 
frequency as our input mode. 



the SQUID ring. In order to verify that the photon oc- 
cupation of the e2,3 modes is due to down conversion 
from mode ei, as stated above, we repeated our calcula- 
tion at bias point A with the system in the initial state 
\ip (0)) = |0) ei ® \a) s ® |0) ea ® |0) ea . The time averaged 
energies (n e2 ) and (n e3 ) for this situation are shown in 



E{A) = S T -S A 



(8) 



This measure has the important property that if it is 
negative definite then A is quantum correlated (i.e. en- 
tangled) with at least some of the rest of the system. 

For this particular example we are interested in the cor- 
relations that exist between the two output mode, e2 and 
e$. The computed entropy, /(e2,es) = S e2 +S e3 — S e2i e 3 , 
is shown in figurc|3fa) as a function of dimensionless time 
uj s t for the bias points A and B in figure Bib). In addi- 
tion, in figure |SJb) we show the individual, and identical, 
entropies of the two (receiver) modes e2 and e3 with the 
initial em field mode e\, i.e. / (ei, e,) = S ei + S ei -S ei;El . 
i = 2 or 3 for the same two bias points A and B. Fig- 
ure EJc) shows a repeat of the calculation of figure Efa) 
having included an additional output mode at the same 
frequency as the input mode, as in figure^fc). Again, the 
inclusion of this extra output mode has modified the re- 
sults. However, the two low frequency output modes (ei 
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FIG. 5: (a) the entanglement entropy as function of dimen- 
sionless time between the two electromagnetic output modes 
of the system of figure |21 calculated at the bias flux points A 
and B in figure |3Jb). (b) the entanglement entropy between 
the input mode and one of the output modes as a function 
of dimensionless time, again calculated at the two bias points 
A and B in figure |3] (c) a repeat of the calculation of (a) 
modified by the inclusion of a third output mode at the same 
frequency as our input mode. 



and 63) are still significantly entangled with each other 
at the appropriate point in external flux, again indicating 
that the excitation of additional output modes reduces, 
but does not suppress, the entanglement in between the 
two low frequency output modes (e 2 and 63). We note 
that the use of an actual LGi£*i£i circuit would altogether 
remove the need to consider higher order harmonics. 

In the photon number expectation values of figure 01 
it is apparent that the decrease in number expectation 
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FIG. 6: (a) the Adami-Cerf entanglement entropies for the 
coupled system of figure|5|as a function of dimensionless time 
calculated at point A in figure|21 where the SQUID ring results 
are in red, the input electromagnetic mode results are in green 
and those for the two output modes are in blue (b) the same 
calculation as above in (a) but now at the external bias flux 
point B in figure El 



value of field mode e\ is directly linked to its increase 
in modes e 2 and e% and that the degree of transfer is 
highly dependent on the value of <& x in the exchange re- 
gion. Correspondingly, from figure EJa) it is clear that, 
with the SQUID ring biased into an exchange region, the 
correlations between the e 2 and e3 modes are greatly af- 
fected by the choice of <1> X . Thus, at B the correlations 
are relatively close to zero at all times while at A it has a 
maximum of around 1.1. This ^-dependent variation in 
the strength of correlation is also seen in figure Efb) be- 
tween the ei and e 2 (or equivalently 63) receiver mode. 
Again these correlations between each pair of coupled 
modes is identical, whether we choose bias point A or B. 
In order to characterise the quantum correlations with 
this system we show in figure [5] the Adami-Cerf entan- 
glement entropies for each mode of the system where, in 
figurcinta), for bias point A in figurcEtb), all four modes 
are entangled with the rest of the system throughout the 
evolution while in figure EJb), while at bias point B, it 
is clear that the level of entanglement is much less. It 
is apparent in these figures that both the von Neumann 
and Adami-Cerf approaches indicate that the level of en- 
tanglement of the modes in the total system can be con- 
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trolled by adjusting the external bias flux on the SQUID 
ring. The control over the entanglement which is man- 
ifest in figure HH is to be compared to the case in quan- 
tum optics. Here, the optical media used are typically 
weakly (polynomially) non-linear as opposed to the ex- 
tremely strong non-linear (non-perturbative) properties 
of the SQUID ring in an exchange region. Furthermore, 
these media generally couple weakly to cm fields which 
does not have to be the case for SQUID rings. Most im- 
portantly, the strength of the ring-field interaction (with 
the resulting two photon down conversion and entangle- 
ment) in an exchange region can be varied by adjusting 
the bias flux on the ring. To the best of our knowledge, 
such control is very difficult to achieve with optical me- 
dia. 



% = \K P\ + E 1 (Mac + 1) (2a sc pal ~ 



where M sc is related to the temperature T = 4.2K, 
and frequency u>t, sc , of each decohering bath for 
each subcomponent (sc) of the system via M sc = 
(cxp (HoJbsc/kBT) — 1) .In equation @ J sc is the cou- 
pling (i.e. the damping rate) between each of the system 
subcomponents and its respective bath. The choice to 
couple each subcomponent to a separate decohering bath 
is in keeping with the standard decoherence model used 
for open quantum systems^. There is, of course, the 
possibility that in certain situations all the elements of 
this system could couple to a mutual decohering bath2&. 
This is an interesting scenario that would require a de- 
tailed analysis beyond the scope of this paper. 

Following on from the entanglement entropy versus 
dimcnsionlcss time plots for the four mode system in 
the absence of dissipation (figure |SJ, we show in fig- 
ures [7] and the subcomponent (Adami-Cerf) entropies 
for two levels of dissipation, namely 7 = 10~ 5 w s and 
7 = 10~ 4 cj s . As can be seen, in the less dissipative 
case shown in figure (which would appear to be at- 
tainable experimentally^^) entanglement is maintained 
over a relatively long time period. However, as one would 
expect, over long enough times any original entanglement 
disappears. Nevertheless, for the case of figure [7J it ap- 
pears that a useful level of entanglement is maintained 
amongst the output modes over a range of over 5000 in 
normalised time or 3 x 10~ 8 s in real time. This is con- 
siderably longer than the time constant corresponding 
to the A S C S oscillator frequency of the SQUID ring, i.e. 
(1/27TVA7CT)' 1 = 7 x 10~ 12 s for A s = 3 x 10~ 10 H and 
C s = 5 x 10~ 15 F. Thus, provided the level of dissipation 
can be kept at or below this level it should prove possi- 
ble to utilise this two photon entanglement in practical 
situations. 



2. With dissipation 

As has been discussed at great length in the litera- 
ture^ coupling of a quantum object to classical dissipa- 
tive environments leads to decoherence over some char- 
acteristic time period. Clearly this is of great importance 
for any realistic discussion of entanglement in multipar- 
tite quantum systems. In order to model the effects of 
dissipation we have adopted the standard approach, very 
often used in the field of quantum optics. In this model 
the components of the system are coupled to decohering 
thermal baths. The master equation for the evolution of 
the density of the multi-component systems then has the 
formS 



I 

III. SIX MODE (FOUR PHOTON) DOWN 
CONVERSION 

The scheme for a six mode system, consisting of one in- 
put oscillator mode (frequency w ei /27r) and four identical 
output modes, all oscillating at one quarter of the input 
frequency, is shown diagrammatically in figure [5] Using 
the same procedure as we adopted in the case of the two 
photon down conversion described above, we have com- 
puted the number expectation values against normalised 
time for the input, the SQUID ring and the four out- 
put modes. These have been computed at the bias flux 
point A' in the centre of the exchange region (divided 
crossing) of figure mil These expectation values are plot- 
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FIG. 7: Repeat of the calculation in figureUJa) but now with 
dissipation introduced for a 7 = 10~°u> 3 with the colour cod- 
ing green (input mode), blue (output mode) and red (SQUID 
ring). 



a\ c a sc p - pa\ c a sc ) + ~^ M sc (2a\ c pa sc - a sc a\ c p - pa sc a\ c ) (9) 
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ted in figure ITTI for the input mode and the four output 
modes, each of which generates identical results. As can 
be seen, in the absence of dissipation the sum of the num- 
ber expectation values for these output modes is reason- 
ably close to 4, indicating that here the down conversion 
process via the SQUID ring is quite efficient. We note 
that the level of efficiency of this down conversion pro- 
cess decreases as the splitting frequency at the divided 
crossing increases. We suspect that this could be cor- 
rected if a sufficiently thorough review of the external 
flux dependence were to be undertaken. However, this 
is impractical to perform given the limitations on our 
current computational power. 

The computational power required to generate accu- 
rate solutions for this four photon down conversion, and 
gain a detailed understanding of the correlations that ex- 
ist between the components of the system, is extremely 
demanding. As with the two photon down conversion 
case, but even more so, the computational power at our 
disposal did not allow us to study the correlations be- 
tween the four (identical) output modes of the system. 
What we have been able to do is to take any of the output 
photons, at a quarter of the frequency of the input pho- 
ton, and show that each is strongly entangled with the 
rest of the system. By appealing to the earlier results pre- 
sented in this paper, it is reasonable to assume that the 
four output modes are entangled with each other in ad- 
dition to any entanglement with the input mode and the 
SQUID ring. The problem of describing, and quantify- 
ing, entanglement for systems of more than two quantum 
particles, is still a matter for serious debat o 38 ! 39 . How- 
ever, in a qualitative way, the computed results displayed 
in figure El suggest a practical route to the entanglement 
of many identical photon states through the intermedi- 
ary of a highly non-perturbative SQUID ring medium. 
Again, with regard to the computational limits we faced, 
we were not able to introduce dissipation into the many 
moded system through a master equation of the type JHJ). 
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FIG. 8: Repeat of the calculation in figureUJa) but now with 
dissipation introduced for a 7 = W~ 4 u> s with the colour cod- 
ing green (input mode), blue (output mode) and red (SQUID 
ring). 




FIG. 9: Schematic for a six mode system with a mesoscopic 
quantum regime SQUID ring inductively coupled to one input 
and four output electromagnetic modes. 



IV. HIGHER ORDER DOWN CONVERSION 
PROCESSES 

Due to the limitations imposed by our computing 
power we have not been able to extend our investiga- 
tions of entanglement beyond the six mode (four pho- 
ton down conversion) range. Nevertheless, following the 
two and four photon processes we now consider a factor 
twenty down conversion process between an input and an 
output em mode. For this calculation we have used the 
same SQUID ring parameters as before. However, with 
this calculation describing a factor twenty frequency con- 
version, the exchange regions have been redistributed as 
shown in figure IT51 For the calculation presented in this 
section we have chosen to flux bias the ring at the cen- 
tre of a particular exchange region in figure fiTA denoted 
by A" . At this value of external flux we have been able 
to demonstrate quantum down conversion by a factor of 
twenty between the input and an output em oscillator 
mode. In order to reduce the computational complexity 
of this calculation, and in contrast to the previous ex- 
amples, we have taken u ei — uj s and assumed an initial 
(start) number occupancy in this input mode of |2) . In 
the computation we have inductively coupled the SQUID 
ring to an output oscillator mode at one twentieth of the 
input photon frequency. Setting the input and output 
coupling strengths at [i = 0.01, we should see, for a per- 
fectly efficient down conversion process, that the number 
expectation value of the output modes reaches a peak of 
approximately 20 as a function of normalised time. The 
actual level of occupancy in the output mode is much 
lower than this potential maximum. However, as dis- 
cussed earlier, the efficiency of this down conversion pro- 
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FIG. 10: (a) the time averaged energies at T = OK for the 
electromagnetic input and output modes of the coupled sys- 
tem of figure [5] as a function of external bias flux with the 
circuit parameters given in the text (b) a zoomed in enlarge- 
ment of the central section of the results in (a) as a function 
of bias flux showing a series of energy exchange regions; in 
the calculations presented subsequently the bias point A' has 
been selected. 
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FIG. 11: Number expectation values for the electromagnetic 
input and output modes of the coupled system of figure [$] as 
a function of dimensionless time, at external flux point A' in 
figure [TU1 starting with a photon occupancy of 1 in the input 
mode. 



cess is greatly affected by the value of the external bias 
flux chosen. The results presented in ngurc lHl arc. for the 
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FIG. 12: The Adami-Cerf entanglement entropies of the for 
the coupled system of figure as a function of dimensionless 
time calculated at point A' in figure 1101 where the SQUID 
ring results are in red, the input electromagnetic mode results 
are in green and those for the two output modes are in blue. 
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FIG. 13: The T = OK time averaged energy for the output 
electromagnetic mode of a factor 20 SQUID ring mediated 
down conversion process as a function of external bias flux 
where the bias point A" in an energy exchange region is used 
in subsequent calculations with the circuit parameters given 
in the text. 



moment, the best that could be obtained within the prac- 
tical time and computational constraints imposed on us. 
In figure ITU the number expectation value for the output 
oscillator mode (e2) is plotted as a function of normalised 
time. This confirms that, at least from a theoretical view- 
point, high order down conversion (and, by implication, 
up conversion) processes can occur between input and 
output cm modes through the intermediary of a quan- 
tum regime SQUID ring. It should be noted that while 
the occupancy of the lower frequency output mode here 
is not very efficient, it is still a non-trivial deviation from 
zero. It is therefore still valid to look upon this result as 
evidence that SQUID rings in the quantum regime can 
mediate high order down conversion, even if, for the mo- 
ment, at low rates of efficiency. By extension, if the com- 
putational power were available to us, it seems reasonable 
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FIG. 14: Number expectation values for the electromagnetic 
input and output modes of the coupled system of figure 1131 
as a function of dimensionless time, at external flux point A" 
in this figure, starting with a photon occupancy of 2 in the 
input mode; also shown in grey is the output mode number 
expectation values for the zero photon occupancy in the input 
mode. 

to assume that (i) we could obtain the same result if we 
were to couple twenty output modes to the SQUID ring 
and (ii) that there is every possibility that these out- 
put modes would be able to entangle with each other. 
Furthermore, given the extremely non-perturbative na- 
ture of a quantum mechanical SQUID ring, we might 
expect that, with the appropriate experimental arrange- 
ments in place, down conversion frequency ratios much 
higher than 20 might be observed. 

V. CONCLUSIONS 

The flexibility afforded by this ^-dependent control 
(tunability) of the maximum coupling strength, and the 
temporal form of the entanglements within ring-field 
modc(s) system, raises interesting possibilities. This con- 
trol might be of great advantage in the preparation of 
qubits in quantum computing and quantum information 
processing^. Our calculations point to a method for 
generating a wide range of levels of entanglement be- 
tween photon states (or, in general, any quantum cir- 
cuits coupled to the SQUID ring) with only very small 
adjustments in bias flux required. This we consider to 
be the important outcome of this work. To give an 
example, for bias point A in the data of figures |5f a) 
and (b) the time evolution shows instances, such as at 
Lu s t pa 10,000, when the entanglement entropy / (e2,e3) 
is high and the entanglement entropy 7(ei,e,), i = 2 
or 3, is low. Even more significantly, in figure we see 



that at dimensionless times such as u> a t ~ 5000 the out- 
put modes are the only modes which are unequivocally 
entangled. This would appear to be of significance in 
many proposed device applicationsi*2*22iiiI where a way 
is needed to use a known input to prepare a set of entan- 
gled photons and then to disconnect the input, without 
significantly disrupting the entanglement so produced. 
We note that in many experimental situationsSLSLSiS 
any readjustment of output entanglement requires the 
development of a new system setup. Many of the ideas 
explored in this work were inspired by concepts in quan- 
tum optics. In fact, there are strong similarities between 
a quantum regime SQUID ring and a device frequently 
employed in quantum optics - the non-linear optical cou- 
pler - this being formed from two or more linearly coupled 
waveguides where at least one is composed of an optically 
non-linear mediuir4L22i£M4i2£. However, there are two 
important differences between optical couplers and the 
SQUID ring: (i) optical non-linear interactions usually 
have a non-resonant character and (ii) optical non-linear 
couplers are difficult to make because of the strict phase 
and frequency matching conditions that must be obeyed 
by all involved processes simultaneously. Hence, it may 
well be that SQUID ring coupled devices, such as those 
presented here, could lead to a feasible way of realising 
(or at least simulating) non-linear optical couplers. 

The manipulation of the different entanglements in the 
model systems considered in this paper, based on just one 
SQUID ring, appears to open up the possibility of cre- 
ating more sophisticated, multi-SQUID ring based, cir- 
cuit networks where the constituent elements could be 
entangled and disentangled as required simply by chang- 
ing bias fluxes. Obviously further work is called for to 
determine how finely entanglement properties can be ma- 
nipulated in more complex systems. As a pointer to the 
use of such complex systems, we have already shown in 
this work that flux controlled, two photon down conver- 
sion {uj\ — > a->2,3 where oj\ = 2^2,3) and entanglement 
can occur in a SQUID ring mediated system as can four 
photon down conversion (ei — > £2.3,4.5). Given the cosine 
nonlincarity of the ring, we would expect the entangle- 
ment of even more photon states to be possible although 
to demonstrate this is beyond the limits of the current 
computational power available to us. 
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